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ABSTRACT 

We present the first statistical analysis of 27 Ultra-violet Optical Telescope (UVOT) 
optical/ultra- violet lightcurves of Gamma- Ray Burst (GRB) afterglows. We have 
found, through analysis of the lightcurves in the observer's frame, that a significant 
fraction rise in the first 500s after the GRB trigger, that all lightcurves decay after 
500s, typically as a power-law with a relatively narrow distribution of decay indices, 
and that the brightest optical afterglows tend to decay the quickest. We find that the 
rise could either be produced physically by the start of the forward shock, when the jet 
begins to plough into the external medium, or geometrically where an off-axis observer 
sees a rising lightcurve as an increasing amount of emission enters the observers line 
of sight, which occurs as the jet slows. We find that at 99.8% confidence, there is a 
correlation, in the observed frame, between the apparent magnitude of the lightcurves 
at 400s and the rate of decay after 500s. However, in the rest frame a Spearman Rank 
test shows only a weak correlation of low statistical significance between luminosity 
and decay rate. A correlation should be expected if the afterglows were produced by 
off-axis jets, suggesting that the jet is viewed from within the half-opening angle Q or 
within a core of uniform energy density Qc- We also produced logarithmic luminosity 
distributions for three rest frame epochs. We find no evidence for bimodality in any 
of the distributions. Finally, we compare our sample of UVOT lightcurves with the 
X-ray Telescope (XRT) lightcurve canonical model. The range in decay indices seen 
in UVOT lightcurves at any epoch is most similar to the range in decay of the shallow 
decay segment of the XRT canonical model. However, in the XRT canonical model 
there is no indication of the rising behaviour observed in the UVOT lightcurves. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) release between 10^^ and 10^* 
ergs during the prompt emission, which lasts from a few 
milliseconds to a few thousand seconds, and is followed by 
an afterglow, which is observed in the X-ray to radio range 
from as little as a few tens of seconds up to several months 
after the GRB trigger. 

The energy is tra nsported in a relativist ic outflow 
l|Meszaros fc Ree3ll997l l that is likely anisotropic l|Sari et al.l 



Il999h and the energy is expected to be released by inter- 
nal a nd external shocks. Internal shocks (|Rees fc Meszaro3 
1 1994 ) are thought to produce the prompt gamma-ray emis- 
sion, while external shocks l|Rees fc Meszaro3 1 19921 ) are 
thought to produce the afterglow. Internal shocks occur 
when shells of material, which are thrown violently from 
the progenitor at different Lorentz factors, overtake each 
other. The external shocks are produced when the shells of 
material are decelerated by the external medium. 

The short duration of the gamma-ray emission and 
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the rapid decay of the afterglow motivated the con- 
struction and launch of Swift, a rapid response satellite. 
Swift houses 3 instruments: the Burst Alert Telescope 
fBAT:lBarthelniv et af] (120051 )'). the X-ray Telescope (XRT; 
iBurrows et al.l ll2005h) and the Ul tra-violet Optical Tele- 
scope (UVOT: lRoming et all l|2005l )'). The energy ranges of 
the BAT and the XRT instruments are 15 keV - 350 keV and 
0.2 - 10 keV, respectively, and the wavelength range of the 
UVOT is I6OOA-8OOOA. The large field of view of the BAT 
(2 str), enables 1/6*'' of the sky to be searched for GRBs at 
any one time. Once a GRB has been detected by the BAT, 
Swift rapidly slews allowing the XRT and UVOT to observe 
the afterglow within a few tens of seconds after the BAT 
trigger. 

Since launch. Swift has produced a large sample of 
UV/optical and X-ray lightcurves which begin soon af- 
ter t he trigger. The high detection rate with the XRT 
(96%; iBurrows etmbOOSl ) allowed a large number of XRT 
lightcurves to be obtained within the first year. The sys- 
tematic reduction of this sample resulted in the discov- 
ery o f a 4 segment canon ical XRT lightcurve (jZhang et al.l 
I2OO6I : iNousek et al]|2006l '). After 2 years of operation the 
UVO T, with a much low er detection rate than the XRT 
(26%; iRoming et al.|[2008h . has detected more than 50 op- 
tical afterglows. This allows for the first time a systematic 
reduction and analysis of a significant sample of GRBs with 
optical afterglows observed with the UVOT and allows an 
investigation of their generic characteristics. 

In this paper we present and analyze a sample of 27 
UVOT lightcurves of GRB afterglows. In Section 2 we ex- 
plain how we selected the sample of UVOT lightcurves 
and in Section 3 we describe how we systematically re- 
duced and analyzed them. In Section 4 we present the re- 
sults and in Section 5 we discuss our findings. Through- 
out the paper we will use the following fiux convention, 
F (X t" v^'^^ with a and /3 being the temporal and pho- 
ton indices respectively. We assume the Hubble parameter 
H(t = 70kms~'Mpc~^ and density parameters f2A=0.7 
and Q.m = 0.3. Unless stated otherwise, all uncertainties are 
quoted at la. 



2 THE SAMPLE 

To investigate the nature of GRB optical/UV lightcurves 
a large number of well sampled, good quality UVOT 
lightcurves were required. The sample was selected ac- 
cording to the following specific criteria: the optical/UV 
lightcurves must be observed in the v filter of the UVOT 
with a magnitude of $;17.8, UVOT observations must have 
commenced within the first 400s after the BAT trigger and 
the afterglow must have been observed until at least lO^s af- 
ter the trigger. These selection criteria ensure the lightcurves 
have adequate signal to noise and cover both early and late 
time evolution. In addition to the above criteria, the colour 
of the afterglows must not evolve significantly with time, 
meaning that at no stage should the lightcurve from a sin- 
gle filter significantly deviate from any of the other filter 
lightcurves when normalized to the v filter. This ensures that 
a single lightcurve can be constructed from the UVOT multi- 
filter observations. Three GRBs, GRB 060218, GRB 060614 



and GRB 060729 were excluded as they showed significant 
colour evolution. 

In total 27 GRBs, which occurred between 1st January 
2005 and 1st August 2007, fit the selection criteria. As there 
were no short GRBs that met the selection criteria, all the 
GRBs in this sample are long. Observations, for the majority 
of GRBs in this sample, began within the first 100s and the 
optical afterglow was detected until at least 10^ s. Formally, 
GRB 050820a meets our selection criteria, but we have ex- 
cluded this burst from the sample because the BAT triggered 
on a precursor of this GRB, and the main GRB took place as 
Swift entered the South Atlantic Anomaly (SAA) with the 
consequence that UVOT completely missed the early phase 
of the afterglow. 



3 DATA REDUCTION 

After the BAT has triggered on a GRB and Swift has slewed, 
the UVOT performs a sequence of pre-programmed expo- 
sures of varying length in multiple observing modes and fil- 
ters designed to balance good time resolution and spectral 
coverage. Observations are performed in either event mode, 
where the arrival time and position are recorded for every 
photon detected, or in image mode, in which the data are 
recorded as an image accumulated over a fixed period of 
time. The pre-programmed observations begin with the set- 
tling and finding chart exposures. The settling exposure is 
not included in this analysis because the cathode voltage 
may still increase during the first few seconds. Two finding 
charts follow immediately after the settling exposure and 
these are observed in event mode with the v and white fil- 
ters (as of the 7th November 2008 the finding charts are 
observed in u and white). The rest of the pre-programmed 
observations are a combination of event and image mode ob- 
servations until ~2700s after the trigger, after which, only 
image mode is used. These observations are taken as a se- 
ries of short, followed by medium and then long exposures, 
which are usually observed with all seven filters of the UVOT 
(white, V, b, u, uvwl, uvm2, uvw2). However, for some tar- 
gets, the pre-programmed observation sequence may change 
or not be executed fully because of an observing constraint. 

For GRB observations, the earliest part of the lightcurve 
is expected to show variability over the shortest timescales. 
As the finding charts, which are exposures of 100-400s, con- 
tain the earliest observations of the GRB, it is essential to 
obtain lightcurves from these event lists as well as obtaining 
lightcurves from the images. 



3.1 Event and Image Data Reduction 

To obtain the best possible lightcurves, we refined the as- 
trometry of the event files before the count rates were ob- 
tained. The astrometry of the event list was refined by ex- 
tracting an image every 10s and cross-correlating the stars 
in the image with those found in the USNO-Bl catalogue. 
The differences in RA and DEC between the stars in the im- 
age and the catalogue were converted into pixels and then 
applied to the position of every event in the event list during 
that particular 10s interval. This process was repeated until 
the end of each event list. 

The images used in this paper were reprocessed by the 
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the Swift Science Data Center (SDC) for the UVOT GRB 
catalogue l|Roming et al.|[2008l ). These images were used be- 
cause not all the image files in the Swift archive have been 
corrected for modulo-8 fixed pattern noise^ ; only those pro- 
cessed with Swift processing script version 3.9.9 and later 
have had this correction applied. For a small number of im- 
ages, the aspect correction failed during the SDC processing 
and in these cases, the images were corrected using in-house 
aspect correction software. 

An aperture of 5", select ed in order to be compatible 
with the UVOT calibration (jPoole et al.ll2008h . was used 
to obtain the source count rates. However, for sources with 
a lo w count rate, it is more precise to use a smaller aper- 
ture (jPoole et al.ll2008l l. Therefore, below a threshold of 0.5 
counts per second the source count rates were determined 
from 3" radius apertures and the count rates were cor- 
rected to 5" using a table of aperture correction factors con- 
tained within the calibration. The background counts were 
obtained from circular regions with radii, typically, of 20". 
These regions were positioned, in each GRB field, over a 
blank area of sky. For each GRB, the same source and back- 
ground regions were used to determine the count rates from 
the event lists and the images. The software used can be 
found in the software release, headas 6.3.2 and the version 
20071106 of the UVOT calibration files. 

For each GRB, to maximise the signal to noise of the 
observed optical afterglow, the lightcurves from each filter in 
which the burst was detected were combined into one over- 
all lightcurve. The lightcurves corresponding to the differ- 
ent filters were normalized to that in the v filter. This was 
possible because there was no significant colour evolution 
between the filters, which was one of the selection criteria 
described in Section 2. The normalization was determined 
by fitting a power-law to each of the lightcurves in a given 
time range simultaneously. The power-law indices were con- 
strained to be the same for all the filters and the normal- 
izations were allowed to vary between the filters. The ratios 
of the power-law normalizations were then used to normal- 
ize the count rates in each filter. The time range used to 
normalize the lightcurves was selected on the basis that the 
time range included data from all filters in which the burst 
was detected and (as far as possible) the lightcurves were 
in a power-law decay phase. Once the lightcurves were nor- 
malized, they were binned by taking the weighted average 
of the normalized count rates in time bins of 5T/T = 0.2. 
The overall lightcurves were converted fr om v count rate t o 
v magnitude using the zero point 17.89 (|Poole et al.ll200^ . 
For many GRBs, the trigger time does not represent the 
true start of the gamma-ray emission. Therefore, the start 
of the gamma-ray emission was chosen as the start time 
of the T90 parameter. This parameter corresponds to the 
time in which 90% of the co unts in the 15 keV - 3 50 keV 
band arrive at the detector l|Sakamoto et al.l |2008| ) and is 
determined from the gamma-ray event data for each GRB, 
by the BAT processing script. The results of the process- 
ing are publicly available and are provided for each trigger 
at http://gcn.gsfc.nasa.gov/swift_gnd_ana.html. The differ- 
ence between the trigger time and the start time of the T90 
parameter is typically less than a few seconds. However, in 
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Figure 1. The rest frame relative flux SEDs for the 21 GRBs for 
which a luminosity lightcurve was produced. The relative flux for 
each filter of each GRB was determined by multiplying the rela- 
tive count rate fro m Table [T] by the flux conversion factor given in 
iPoole et al. I 1I2OO8I '). The relative flux values have been corrected 
for Galactic extinction. The vertical line at I6OOA marks the rest 
frame wavelength at which the luminosity lightcurves were pro- 
duced. 



a minority of cases the difference is much larger, with the 
largest difference being 133.1s for GRB 050730. 

3.2 Photometric Redshifts 

Spectroscopic redshifts were obtained from the literature 
for 19 of the GRBs in the sample (see Table [1]). For a 
further 4 GRBs, it was possible to determine the redshift 
using an instantaneous UVOT-XRT Spectral Energy Dis- 
tribution (SE D), which was created using the method of 
ISchadv et aD ([2008). The SEDs were fit with the best fit- 
ting model of either a power-law or broken power-law, with 
Galactic and host galaxy absorption and extinction. The 
Galactic Nh was taken from t he Leiden/Arge ntinc/Bonn 
(LAB) Survey of Galactic HI (|Kalberla et al.l i2005. ) and 
the Galactic extinction was taken from a composite 
lOO^tr n map of CO BE/DIRBE and IRAS/ISSA observa- 
tions jSchlegel et al . 1998). For the host extin ction the Small 
Magellanic Cloud (SMC) extinction curve jPeil Il992l ) was 
assumed. The host reddening and absorption, E(B-V), Nh, 
and the redshift were left to vary. The resulting redshifts 
can be found in Table |2] 

3.3 Luminosity Lightcurve 

Luminosity lightcurves were produced for all the GRBs 
whose host E(B-V) value could be determined, except for 
GRBs with photometric redshifts which have a la error on 
the redshift that corresponds to an uncertainty in log lumi- 
nosity of > 0.1. In total, luminosity lightcurves were pro- 
duced for 21 of the GRBs in the sample. 

For the 21 GRBs, the observed count rate was converted 
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into luminosity at a common rest frame wavelength. In or- 
der to select the common wavelength and determine the re- 
sulting k-correction factor for each lightcurve, an SED was 
computed for each GRB. The SED was produced by mul- 
tiplying the relative count rates in each filter, given in Ta- 
ble [J by th e coun t rate-to-flux conversion factors given by 
iPoole et al.l (120081 ) . These relative flux densities were cor- 
rected for Galactic extinction and positioned at the central 
wavelength of the filter in the rest frame of the GRB; the 
SEDs are shown in Fig. [1] The common rest frame wave- 
length at which to determine the luminosities was selected to 
maximise the number of GRBs with SEDs that include this 
wavelength and to be relatively unaffected by host extinc- 
tion. The wavelength which best satisfies these conditions is 
I6OOA. The k-correction factor was taken as the flux density 
in the rest frame at I6OOA, -Fieoo, divided by the flux density 
at the observed central wavelength of the v filter (5402A), 
Fv. In the case where a GRB's SED did not cover I6OOA, an 
average k-correction was determined from the other GRBs 
in the sample, whose SEDs covered both I6OOA and the rest 
frame wavelength corresponding to the v filter. 

To produce the luminosities, the lightcurves in count 
rate were corrected for Galactic extinction, converted into 
flux density and then into luminosity using the following 
equation: 

L(1600) = 4TvDlF^k (1) 

where L(1600) is the luminosity at a I6OOA, Dl is the lu- 
minosity distance. A; = (1 -|- 2;)(-Fi60o/fv) is the k-correction 
factor, and z is the redshift of the GRB. Finally, the luminos- 
ity lightcurves were corrected for host extinction using the 
^1600 values given in Table [1] These values we re determined 
for th e GRBs, using the Av values reported in lSchadv et al.l 
l|2008l ). 

3.4 Bolometric Energy and GRB Classification 

The k-corrected isotropic energy of the prompt gamma-ray 
emission E^jso, was calcul ated for each GRB with known 
redshift using Eq. 4 from iBloom et al.1 12001,; '). The ener- 
gies were corrected to a rest frame bandpass of 0.1 keV to 
10000 keV using one of three spectral t emplates: a power - 
law, cut-off power- law or Band function l|Band et al.lll993 ). 
As Konus-Wind has a larger energy range than BAT, the 
spectral analysis of the prompt emission, observed by Konus- 
Wind, better represents the spectral behaviour. These re- 
sults were taken from the literature. The spectra observed 
with Konus-Wind tend to be best fit by a cut-off power- 
law, therefore the spectral template chosen for the GRBs 
that were not observed with Konus-wind, with power-law 
spectra of photon index F > —2 in the 15 keV - 150 keV 
energy range w as a cut-off power-law with Epeak=W2.2 keV 
l|D'Alessio et al. 2006). For the GRBs with a power- law 
spectrum with a photon index of F < —2 in the 15 keV - 
150 keV en ergy range, a Band function was used with 
Fi = -0.99 (|D'Alessio et al.ll2006l ) and Epeak=15 keV. The 
resulting k-corrected energies can be found in Table [3] 

The GRBs in this sample were classified into three cate- 
gories depending on the ratio R of the fluence in the 25 keV - 
50 keV and 50 - 100 keV B AT energy bands, which are given 
in ISakamoto et al.1 ||2008D . The cate gories and their respec- 
tive ratios are: an X-ray flash (XRF) for R > 1.32; an X-ray 



rich GRB (XRR) for 0.72 < 7? ^ 1.32; or a classical GRB 
(C-GRB) for R sC 0.72. Table H lists the GRBs with their 
classifications. In total, there are 12 C-GRBs, 13 XRR and 
1 XRF, which is GRB 060512. 



4 RESULTS 

In this Section, the flux lightcurves in the observer frame 
and their properties shall be investigated first, then the lu- 
minosity lightcurves in the rest frame and their properties 
shall be examined. 

4.1 Observer Frame Flux Lightcurves 

The lightcurves, shown in Fig. [S] are ordered by peak mag- 
nitude from the brightest, GRB 061007, to the faintest, 
GRB 050712. Data points with signal to noise below 2 are 
shown as Scr upper limits. The peak magnitude was taken 
as the maximum magnitude in each binned lightcurve and 
is given in Table |3] A trend is observed in the figure, with 
the brightest GRBs decaying more quickly than the faintest 
GRBs. The lightcurves generally follow one of two types of 
behaviour. Either they rise to a peak within the first 1000s 
and then decay, or they decay from the beginning of the 
observations. There are 6 lightcurves that appear to rise to 
a peak between 200s to 1000s after the burst trigger. The 
peak times for the 6 GRBs were determined from a Gaus- 
sian fit to each lightcurve in log time. The fit was performed 
between the brightest data point and the data points on ei- 
ther side to the point at which the count rate is 60% of the 
peak value. The mean of these peak times is 397s. For the 
remaining GRBs, the beginning of the lightcurve was taken 
as the upper limit to the peak time and the mean of these 
upper limits is 132s. 

To classify the behaviour of the lightcurves, two power- 
laws were fit to each lightcurve, covering the time ranges 
from the start of observations until 500s and from 500s until 
the end of the observations; the results are given in Tabled 
A time of 500s was chosen as it ensures that the early power 
law fits are performed on at least 100s of each lightcurve 
and because the rising phase tends to occur during the first 
500s. 

A comparison of the temporal indices before and af- 
ter 500s is shown in Fig. O The figure is divided into three 
groups of behaviour, which are: lightcurves which rise be- 
fore 500s, lightcurves which decay more steeply after 500s 
and lightcurves which decay less steeply after 500s. Each 
of these groups contains a similar number of GRBs. From 
Fig. [S] there are 4 lightcurves that are clearly rising with an 
Q<5oos ranging from 0.26±0.13 (GRB 060605) to 0.73±0.14 
(GRB 070420). A further 7 lightcurves are consistent with 
Q!<500s ~ or have large errors and thus it is not clear 
if these lightcurves are rising, constant or decaying before 
500s. The remaining 15 lightcurves decay with temporal 
indices of between a^soos = -0.12 ± 0.05 (GRB 061121) 
and a<5oos = -2.67 ± 0.80 (GRB 050726). After 500s, aU 
the lightcurves decay with values of a>5oos ranging from 
-0.50 ±0.05 (GRB 050712) to -1.67 ±0.15 (GRB 070420), 
except for GRB 050726 where, due to the poor signal to 
noise, it is not possible to tell if the lightcurve after 500s is 
rising or decaying. 
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Figure 2. The lightcurves, ordered by peak magnitude from brightest (GRB 061007) to faintest (GRB 050712). With a few exceptions, 
the hghtcurves appear also to be ordered by decay rate with the brighter bursts decaying more rapidly than the fainter bursts. The down 
arrows in each lightcurve represent 3cr upper limits. 



The mean and intrinsic dispersion of the temporal in- 
dices was determined u sing the maximum likelihood method 
l|Maccacaro et allllQSSl ) , which assumes a Gaussian distribu- 
tion. The mean for q:<5oos is — 0.48to'i9 with a dispersion 
of 0.69to'o6 and the mean for q:>5003 is — 0.881q;o7 with a 
dispersion of O.SIIq qj. To see if q:<500s and a>5oos are inde- 
pendent parameters, a Spearman rank test was performed. 
This test gives a coefficient of -0.22 with a probability of 
73%, indicating no evidence for a correlation and suggest- 



ing that the behaviour after 500s is independent of the be- 
haviour before 500s. 

Since the lightcurves in Fig. [2] suggest a connection be- 
tween the brightness and the decay rate, a Spearman rank 
correlation was performed between the temporal indices and 
the interpolated magnitudes of the lightcurves at 400s. The 
test performed between q:<5oos and the magnitude at 400s, 
indicates that these parameters are not related, as the co- 
efficient is -0.28 at 84% confidence. However, the Spearman 
rank test performed on q>5oos and the magnitude at 400s 
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57 


3 


42 





92 










N/A 



Table 1. Spectroscopic redshifts were largely taken from the literature. For four GRBs, photometric redshifts, indicated by an *, were 
determined using the X RT-UVOT SEDs (s ee Section 3 for more details). The host extinction values at I6OOA, were calculated from 
the best fit Av given in llSchadv et al.|[200^ . The count rate relative to v is provided for each filter detected by the UVOT. For a few 
GRBs, there were no observations with the wh ite filter an d for G R B 060804, the whit e d ata could not be u sed ; the white column for 
these GRBs coiitain N/ A. References: a^ljako bsson ct all ll2006al) b) iFoley et al.l ||2005 ) c) IChen et al.l ||2005^ d) Fvnbo et all (120051') e) 
I Jakobsson et al.l ll2006ah f ) iQuimbv et al.l ll2005i) g) iFvnbo et al l 1I2OO6I) h) Bereer et al, | |2006^ i) jProchaska et"aL ll200d) il iBloom et alj 
(I2OO6I') k) jjakobsson et al.l l2006a^ DiPeterson fc Schmidt] 1I2OO6I) m'llLedoux et al.l 1I2OO6I) n'llRol et al.l 1I2OO6I ') o) jjakobsson et al.l 1 I2OO6I ') 
I Jakobsson et"^ 1 1|2006^ q lBloom et al. I ||2006|'| r) I Jaunsenetah I ll2007h sl lBereer et al.1 ||2007|) 



gives a coefficient of 0.59 at 99.8% confidence (see Fig. U}. 
The correlation is statistically significant at > 3cr and there- 
fore implies that brighter GRBs tend to have faster decays. 

The mean redshifts of the two columns in Fig. [2] are 
(z) = 1.63 and {z) — 2.98, suggesting that the decay rate 
and magnitude are also correlated with redshift because the 
optical afterglows that are brighter and decay more steeply 
tend to have lower redshifts. A Spearman Rank test per- 
formed between a>5oos and redshift gives no evidence for 
a correlation with a correlation coefficient of 0.07 and a in- 
significant probability of 23%. However, a Spearman Rank 
test performed between the redshift and magnitude at 400s, 
provides a weak correlation with a coefficient of 0.32 and 
a statistical significance of 87%. Moreover, if the Spear- 
man Rank test is performed between the peak magnitude 
and redshift, the link between the redshift and magnitude is 
stronger and more significant as the coefficient is 0.55 and 
the confidence is 99.3%. These correlations imply that the 
correlation between magnitude at 400s and a>5oos is only 
weakly dependent on redshift. 

The lightcurves of the GRBs after 500s, in a few cases, 
appear to show a change in their temporal behaviour. To 
quantify this behaviour, a broken power-law was fitted to 
each GRB from 500s until the end of the observations. 
The broken power-law is considered an improvement if the 
X^/D.O.F has decreased and the probability of chance im- 
provement is small (< 1%), as determined using an F-test. 
In five cases a broken power-law was an improvement com- 
pared with a single power-law. The results of the broken 
power-law fits for the 5 GRBs are given in Table [5| In four 
of these five cases, the broken power-law shows a transi- 
tion from a shallow to a steeper decay. In the fifth case, 
GRB 070318, the decay became shallower at late times. 
To test if a single break is sufficient for the decay after 



GRB 


Photor 


QGtric Rcdshift 


X^/D.O.F 


GRB050801 


1. 


38 ± 0.07 


(14/12) 


GRB06070S 


1. 


92 ± 0.12 


(21/20) 


GRB061021 





....+0.06 

-0.01 


(229/174) 


GRB070420 


3 


n,+0.96 
-0.68 


(58/60) 



Table 2. Photometric redshifts for 4 of the GRBs without spec- 
troscopic redshifts. 



500s, a doubly broken power-law was fit to these 5 GRBs. 
As with the broken power-law, the doubly broken power- 
law was considered an improvement if the x^/D.O.F de- 
creased and the probability of chance improvement is small 
(< 1%). The doubly broken power-law was an improve- 
ment for only GRB 070318. The best fit values for this 
model are: ai = -1.08 ± 0.01, hreak,! = 538001 ^^EJ^J , 
aa = -O.ntoll, tbreak,2 = 197000l , fts = "1.72±0.18 
with x^/D.O.F = 62/23. 

For the 4 GRBs where the broken power-law was the 
best fit, the mean decay index before the break is —0.60 ± 
0.14 with a dispersion of 0.19lQ|!,y and the mean decay index 
after the break is -1.53±0.19 with a dispersion of 0.22lggg. 
The break times range from 6000li5^gs to (4.9±0.25) x 10*s. 
If the mean decay index after 500s is recalculated including 
only those lightcurves that decay as a single power-law after 
500s, the mean is -0.87to Jg with a dispersion of Q.'i^t^oM- 
This mean is similar to the mean decay index determined 
using all the GRBs in the sample. For the lightcurves that 
show a break, the mean decay before the break is consistent 
within 2(T with these mean values. 
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Figure 3. Temporal decay after 500s versus temporal decay be- 
fore 500s. The dashed (red) line indicates the point at which 
the temporal index before 500s equals the temporal index after 
500s. The dashed dotted (green) line indicates the cut off between 
GRBs that rise within the first 500s (above the line) and GRBs 
that decay within the first 500s (below the line). 



4.2 GRB Rest Frame Luminosity Lightcurves 

The luminosity ligiitcurves at I6OOA, in units of erg s~^ 
cm~^ , before and after correction for liost extinc- 
tion, are siiown in Fig. [5] Panel (a) shows the luminosity 
lightcurves before any correction for the host extinction has 
been applied. In both panels of Fig.[5l GRB 060512 lies sig- 
nificantly below all the other lightcurves. We suspect that 
this is caused by either an incorrect determination of the 
host extinction or of the redshift. The redshift of this GRB 
could be wrong because it was not determined from the af- 
terglow spectra , but was based on the alignment of the GRB 
with a galaxy l|Bloom et al.ll2006l ). In lSchadv erall l|2008l ). 
the best fitting model to the SED of GRB 060512 gives a 
poor fit with a /D.O.F = 84/23 and a host extinction of 
E{B — V)host ~ 0.16^0 0^. A photometric redshift was deter- 
mined for this GRB using the method described in Section 



2, of 2 = 2.279 



+0.09 



and an extinction of Ai 



0.00: 



-0.18 "'^^^ "'^^ cAuiii^unjii wi ^1600 — >J.>JU_Q OQ. 

Using these values a luminosity lightcurve for GRB 060512 
at I6OOA was produced. This photometric redshift changes 
the rest frame relative fiux SED and consequently the k- 



I 



O 



o 



14 16 18 

Magnitude 400s after the trigger 



20 



Figure 4. Magnitude at 400s against the temporal index after 
500s, a>500s. A Spearman rank correlation test performed on 
these parameters gives a coefficient of 0.59 at a statistical signif- 
icance of 99.8%, indicating that these two parameters are corre- 
lated. 



GRB 








( 


cm) 




f^2 






050319 


BAND 


8 


60 




-|^q28 


-0-99 


-2.02" 


162.0 


1.44c-f 53 


050525 


CPL 


1 


10 




j^q28 


-1.10'' 


-2.31 


84.1^ 


1.59C+53 


050730 


CPL 


1 


04 




1029 


-1.53" 


-2.31 


162.2 


3.85e+54 


050801 


CPL 


3 


04 


X 


10^® 


-1.99*^ 


-2.31 


162.2 


7.67e+53 


050802 


CPL 


3 


96 


X 


10^® 


-1.54°' 


-2.31 


162.2 


5.86e+53 


050922c 


CPL 


5 


38 


X 


10^® 


-1.55'' 


-2.31 


162.2 


4.30e+54 


051109a 


CPL 


5 


82 


X 


10^® 


-1.25*^ 


-2.31 


161.0*^ 


6.43e+53 


060206 


CPL 


1 


12 


X 


1029 


-1.20" 


-2.31 


78.0*^ 


3.42e+53 


060223a 


CPL 


1 


24 


X 


1029 


-1.74" 


-2.31 


162.2 


4.47e+54 


060418 


CPL 


3 


34 


X 


1028 


-1.50"= 


-2.31 


230.0^ 


2.94e+54 


060512 


BAND 


7 


56 


X 


1027 


-0.99 


-2.48" 


162.2 


5.04O+50 


060526 


BAND 


8 


54 


X 


]^q28 


-0.99 


-2.01" 


162.2 


1.37c4-53 


060605 


CPL 


1 


04 




1029 


-1.55" 


-2.31 


162.2 


1.15e+54 


060607a 


CPL 


8 


10 


X 


-|^q28 


-1.47" 


-2.31 


162.2 


1.81e+54 


060708 


CPL 


4 


56 




-Lq28 


-1.68" 


-2.31 


162.2 


3.79e+53 


060908 


CPL 


6 


08 


X 


10^® 


-1.00" 


-2.31 


151.0" 


1.88e+53 


060912 


POWER 


1 


88 


X 


10^® 




-2.31 


162.2 


2.52e+54 


061007 


BAND 


2 


72 


X 


10^® 


-0.70f 


-2.619 


399. of 


2.00e+53 


061021 


CPL 


1 


48 


X 


-|^q28 


-1.22^ 


-2.31 


777.0^ 


9.65e+52 


061121 


CPL 


2 


86 


X 


-Lq28 


-1.32^ 


-2.31 


606.0^ 


2.75e+54 


070318 


CPL 


1 


63 


X 


10^® 


-1.42" 


-2.31 


162.2 


8.78e+52 


070420 


CPL 


7 


87 


X 


10^® 


-1.23^ 


-2.31 


147. 0-? 


6.39e+54 


070529 


CPL 


6 


14 


X 


10^^ 


-1.34" 


-2.31 


162.2 


6.46C+53 



Table 3. Properties of the GRBs with spectroscopic or pho- 
tometric redshifts. This table contains the luminosity distance, 
the gamma-ray photon indices and peak energies used to de- 
termine the k-corrected isotopic energy for each GRB in the 
energy ra nge 10 keV - 10 MeV. References: a' Sakamot o et al.l 



poos), b'Golenctskii et al.' ('2005a'), c^Golenetskii et al. (2005 




e)Golenctskii et al. 
Golcnctskii^t_alJ 
Golenetskii et al.l 




d ' Golc nctskii ct al. 
f ]Golenetski^_et_alJ 
h lGolenetskilet al.l 
i lGoIenetskii et al ] (|200"7t ^ 



correction factor. The result is that the luminosity lightcurve 
increases by ~ 3 orders of magnitude, which means that this 
GRB is no longer separated from the rest of the GRBs in the 
sample. Nonetheless, this GRB may be intrinsically different 
to all the other GRBs in the sample as this GRB is the only 
XRF in the sample and it may be that XRFs are a class of 
sub-luminous GRBs. However, as it is uncertain whether the 
redshift of GRB 060512 is correct, the luminosity lightcurve 
for this GRB will be excluded from any further analysis. 
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To produce luminosity distributions, the luminosities 
were interpolated from the lightcurves before and after cor- 
rection for host extinction at the 3 rest frame epochs: 100s, 
1000s and 10 ks. The logarithmic distribution of the lumi- 
nosities at the three epochs are shown in Panels (a) to (f) 
of Fig. [S] The distributions consisting of the luminosities at 
100s contain 18 GRBs whereas the distributions for the lu- 
minosities at 1000s and 10 ks contain 20 GRBs. Panels (a) to 
(c), show the logarithmic distribution of luminosities before 
correction for host extinction. The means of these distribu- 
tions at 100s, 1000s and 10 ks in the rest frame are 11.08, 
10.29 and 9.39, respectively. The standard deviations at the 
three rest frame epochs are 0.65 at 100s, 0.71 at 1000s, and 
0.68 at 10 ks. Panels (d) to (f) show the logarithmic dis- 
tributions of the luminosity lightcurves that have been cor- 
rected for host extinction. The mean of the host extinction 
corrected distributions at 100s, 1000s and 10 ks in the rest 
frame are 11.29, 10.55 and 9.64, respectively and the stan- 
dard deviations are 0.57 at 100s, 0.67 at 1000s, and 0.62 at 
10 ks. 

The distribution of rest frame peak times and upper 
limits is shown in Fig. [7] The peak times of the GRBs with 
observed peaks overlap with the upper limits of the GRBs 
without observed peaks. Therefore, it is not possible to tell 
if the GRBs with observed peaks are a separate class, or 
if they belong to the tail end of a distribution where the 
majority of GRB peaks occur before the UVOT can observe 
them. 

To determine if the relationship between the brightness 
of the afterglow and the late time decay rate is intrinsic, 
the luminosity lightcurves were fitted with a power-law from 
150s onwards Q>i50s,res*, where 150s ~ 500s/(l -I- (z)) and 
(2:) = 2.21 is the mean redshift of the GRBs in the sample, 
and a Spearman rank test was performed between this decay 
and the extinction corrected luminosity at 100s in the rest 
frame. These parameters are shown in Fig. |8] The Spearman 
rank test does not support or refute a correlation between 
these parameters because the coefflcient is -0.29 and the 
probability of correlation is not significant at 76%. 



5 DISCUSSION 

In this section, we shall discuss the possible mechanisms 
that could produce the rising behaviour of the early af- 
terglow viewed in the observer frame before 500s, and we 
shall discuss the late afterglow from 500s onwards. We will 
also discuss the implications of the luminosity distribution 
and compare the UVOT lightcurves with the XRT canonical 
lightcurve model. 



5.1 The early UVOT afterglow 

There are several physical mechanisms and geometric sce- 
narios that may produce a rise in the early optical afterglow. 
In this section the following mechanisms and scenarios shall 
be discussed: (i) passage of the peak synchrotron frequency 
of the forward shock Vmj, through the observing band, (ii) 
a reverse shock, (iii) decreasing extinction with time, (iv) 
the onset of the forward shock in the cases of an isotropic 
outflow or a jet viewed in a region of uniform energy den- 



GRB 


Lorcntz Factor at Peak 






050319 


>168 


<2.40c-04 


<9.33c + 16 


050525 


>179 


<2.49c-04 


<9.25c + 16 


050730 


174 


6. 18c-03 


2.72c + 17 


050801 


>275 


< 7. 81c- 04 


<1. 17c 4- 17 


050802 


> 156 


< 1.05c-03 


<1.56c + 17 


050922c 


>308 


<3.91c-03 


<1.93c + 17 


051109a 


>236 


<7.61c-04 


<1.22c + 17 


060206 


>338 


<2.83c-04 


<7.80c + 16 


060223a 


>409 


<3.06c-03 


<1.62e + 17 


060418 


193 


4.27C-03 


2.32C + 17 


060512 


> 72 


<1.95c-06 


<2.48c + 16 


060526 


>247 


<1.55c-04 


<7.09c + 16 


060605 


177 


1.82C-03 


l.SOc + 17 


060607a 


220 


2.30C-03 


1.81C + 17 


060708 


>256 


<4.14c-04 


<9.72e + 16 


060908 


>231 


<2.27c-04 


<8.23c + 16 


060912 


>235 


<3.00c-03 


<1.94c4-17 


061007 


>126 


<4.42c-04 


<1.26c + 17 


061021 


>173 


<1.56c-04 


<7.99c4-16 


061121 


>337 


<2.28c-03 


<1.57g4-17 


070318 


103 


2. 38c- 04 


1.09C + 17 


070420 


228 


7.84C-03 


2.69C + 17 


070529 


>234 


<7.71c-04 


<1.23c4-17 



Table 6. Properties of the GRBs, derived assuming the rise of 
the forward shock is the cause of the rise observed in the UVOT 
Hghtcurves. The initial Lo rentz factors were dete rmined using the 
peak times and Eq. 1 of iMolinari et al] ||2007|) . Where only an 
upper limit to the peak time is known, only a lower limit to the 
Lorentz factor is given. The fraction of mass as baryons and the 
deceleration radius are determined using the Lorentz factor. 

sity, (v) the rise produced by an off-axis jet, which may be 
structured, and finally, (vi) a two component outflow. 

5.1.1 Passage of the synchrotron frequency, i^mj 

The flrst mechanism, the passage of the peak frequency of 
the forward shock which moves with time as I'm,/ oc t~^^^, 
through the observing band is expected to produce a chro- 
matic peak in the optical lightcurve evolving from the short- 
est wavelengths through to the longest wavelengths. For 5 
of the lightcurves with a peak, the UVOT observed the ma- 
jority of the rise and the peak during the two flnding chart 
exposures observed in white and v (see Fig. |9} . If the peak 
were due to the passage of the synchrotron frequency, a 
stepped decrease in flux would be observed in the normalized 
lightcurves at the transition between the white and v obser- 
vations, which has not been observed in any of these GRBs. 
For the 6th GRB with an observed peak, GRB 050730, the 
rise was observed in the v and b fllters. If the passage of j 
was the cause of this rise, the afterglow would appear to be 
brighter in the b filter than in the v filter during the rise, 
and begin to decay earlier. However, the afterglow is not 
brighter in the b filter than in the v filter during the rise. 
Therefore, the passage of Vm,f through the optical band is 
not responsible for any of the peaks observed in these optical 
afterglows. 

During the passage of Vm from the shortest wavelengths 
through to the longest wavelengths, the spectrum of the op- 
tical afterglow changes (assuming slow co oling ) from ly^^'^, 
for u < Um to for Vm < v < Vc ((Sari eir'alJll998l 'l. 

As Vm passes from high frequencies to lower frequencies, 
there will be a change in colour. The colour change between 
the white and v filters can be calculated using the central 
wavelengths of the white and v filters, given in iPoole et al.l 
(l2008h . converted in to frequency: i^whitc = 8.64 X lO^"* Hz, 
Vv = 5.49 X lO^'^Hz and assuming p = 2.3 where p is the elec- 
tron energy index. The colour change between white and v 
as Vm moves from above the white frequency to below the 
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Name 


Classification 


Peak A-Iagnitudc 


Peak Time (s) 




(x^/-D.O.F) 


s 


(x^/-D.O.F) 


050319 


XRR 


17.09 


< 234 


0.09 ± 0.35 


0/2 


-0.63 zb 0.03 


34/17 


050525 


XRR 


13.57 


< 78 


-1,25 ± 0.03 


21/ 9 


-1.00 zb 0.01 


398/26 


050712 


C-GRB 


17.77 


< 178 


0.11 ± 0.64 


4/ 3 


-0.50 ± 0.05 


52/26 


050726 


C-GRB 


17.21 


< 159 


-2.67 ± 0.80 


1/ 3 


0.13 zb 0.29 


21/13 


050730 


C-GRB 


17.22 


744 


0.15 zb 0.50 


6/3 


-0.89 zb 0.05 


72/16 


050801 


XRR 


15.26 


< 66 


-0.46 ± 0.04 


43/10 


-1.17 zb 0.03 


41/15 


050802 


XRR 


17.07 


< 289 


0.07 zb 0.48 


0/1 


-0.75 zb 0.01 


53/29 


050922c 


C-GRB 


14.34 


< 109 


-1.01 ± 0.05 


21/ 6 


-0.94 zb 0.01 


88/14 


051109a 


C-GRB 


16.33 


< 122 


-0.47 ± 0.45 


0/ 2 


-0.68 zb 0.02 


80/16 


060206 


XRR 


16.64 


< 57 


-0.02 ± 0.50 


9/ 4 


-0.70 zb 0.03 


205/22 


060223a 


XRR 


17.33 


< 88 


-1.06 ± 0.33 


4/ 6 


-1.53 zb 0.60 


31/14 


060418 


XRR 


14.69 


260 


0.01 ± 0.03 


138/ 6 


-1.22 zb 0.01 


58/21 


060512 


XRF 


16.50 


< 114 


-0.74 ± 0.08 


3/ 6 


-0.99 zb 0.05 


10/14 


060526 


XRR 


16.59 


< 82 


-0.29 ± 0.08 


36/ 8 


-0.69 ± 0.02 


141/22 


060605 


XRR 


16.50 


459 


0.26 ± 0.13 


6/ 4 


-0.94 zb 0.03 


27/11 


060607a 


C-GRB 


14.50 


254 


0.50 zh 0.03 


702/ 8 


-1.17 ± 0.03 


152/20 


060708 


XRR 


17.19 


< 72 


-0.10 zb 0.10 


8/ 9 


-0.89 ± 0.03 


39/23 


060804 


XRR 


17.43 


< 231 


-1.76 zb 0.49 


0/2 


-0.57 ± 0.08 


21/9 


060908 


C-GRB 


15.21 


< 88 


-1.19 zb 0.05 


6/ 8 


-1.03 zb 0.04 


20/18 


060912 


XRR 


16.44 


< 114 


-0.97 zb 0.09 


2/ 6 


-0.53 zb 0.01 


149/27 


061007 


C-GRB 


12.68 


< 298 


-1.63 ± 0.11 


1/ 1 


-1.59 zb 0.01 


129/28 


061021 


C-GRB 


15.64 


< 79 


-0.92 zb 0.06 


3/ 9 


-0.88 zb 0.01 


556/29 


061121 


C-GRB 


15.67 


< 53 


-0.12 zb 0.05 


181/11 


-0.67 zb 0.01 


179/31 


070318 


C-GRB 


15.37 


316 


0.42 zb 0.03 


19/ 9 


-1.02 zb 0.01 


261/27 


070420 


C-GRB 


17.16 


347 


0.73 zb 0.14 


13/ 5 


-1.67 zb 0.15 


56/23 


070529 


C-GRB 


15.95 


< 131 


-1.64 zb 0.14 


5/ 5 


-0.53 zb 0.04 


42/19 



Table 4. The v peak magnitude, peak time and best fit decay indices of the lightcurves before and after 500s and the classification 
of each GRB. In the cases where a peak was not observed an upper limit to the peak time is provided. The GRBs are classified into 
three categories depen ding on the ratio of the fleuence in the 25 keV - 50 keV and 50 - 100 keV BAT energy bands, which are given in 
ISakamoto et al.l ||2008|) : an X-ray flash (XRF), an X-ray rich GRB (XRR) or a classical GRB (C-GRB). 



GRB 



(X^/D.O.F) Ax^ 



GRB 050525 
GRB 050922c 
GRB 060526 
GRB 061021 
GRB 070318 



-0.80 ± 0.01 
-0.76 ± 0.03 
-0.33 ± 0.04 
-0.51 ± 0.02 
-1.09 ± 0.01 



6000^ 



,,,ool-» 

n+1100 
-1000 
30800^-00 

49300+2500 



-1.70 ± 0.08 
-1.20 ± 0.05 
_2 33+0.65 
-0.47 

-0.78 ± 0.03 



( 120/24 ) 
( 16/12 ) 
( 32/20 ) 
( 49/27 ) 

( 131/25 ) 



93 
507 
129 



Table 5. Five of the lightcurves in the sample could be fit with a broken power-law from 500s onwards. The final column provides the 
difference in (for 2 additional degrees of freedom) when fitting the lightcurve after 500s with a single and broken power laws. 



V frequency is 0.48 magnitudes. In the lightcurves of Fig. |9l 
which have been normalized using the late time data, the 
colour difference between the white and the v filter after the 
peak is zero. None of the 5 lightcurves, in Fig. [O] with early 
white and v observations, show such a large offset between 
the lightcurves in two filters during the rise. Furthermore, 
th e equation for Vmj i n a constant density medium, as given 
bv lZhang et"al] l|2006l ) as; 

Vm.f = (6 X 10'"H^)(1 + zf'^El!,^eyji''{t/lday)-'"^ (2) 

where E52 = 10^^ E is the isotropic energy in units of 10^^ 
ergs, Ee is the fraction of energy in the electrons, es is the 
fraction of energy in the magnetic field, z is the redshift and 
t is the time. Assuming ee and es are not evolving with time 
and given the time t\ at which Umj is at a given frequency 
v\, the time t2 at which Umj is at the frequency 1/2 is given 
by: t2 = t\(y2li'\f~'^^'^'' ■ Using the central wavelengths of 
the white and v filters, the v filter should peak 1.35x later 
than the white filter. There does not appear to be any time 
difference between the peak in the white and v filters in the 5 
Ughtcurves shown in [51 therefore, the passage of Vra through 
the optical band is not the cause the rise in the optical band. 

For I'm,/, to have passed below the v filter (5402A) by 
the time the UVOT has begun observations it ~ 100s), using 
(£^fe,iso) = 1.5 X lO^'* erg and (z) = 2.2 1, then fjej^'^ < 1.7 x 
10~^. The values ee and es provided in lPanaitescu fc Kumaij 
(|2002l ) for 10 GRBs give values for €^eg ranging from 3 x 



10 ^ to 2 x 10 , suggesting that 1.7 x 10 ^ is consistent 
with values found for other GRBs. 



5.1.2 The reverse shock 

Considering a constant density medium, there are two main 
cases of the reverse shock that depend on the position of 
the peak synchrotron frequency of the reverse shock i^m,r 
relative to the optical band Vopt. If i^m,r < Vopt then the 
lightcurve produced by the reverse shock is expected to 
decay immediately after the peak with a = (3p + l)/4 
(|Zhang et al.ll2003l ). where p is the electron energy index. 
The value of p is typically taken to lie between 2 and 
3, therefore a is expected to range between a = —1.75 
for p = 2 to a = —2.5 for p = 3. Within the sample, 
GRB 050726, GRB 061007 and GRB 070529 are the only 
GRBs with decays before 500s that are consistent at 2cr with 
the slowest expected reverse shock decay index a — —1.75; 
GRB 050726 has a decay of a = -2.67 ± 0.80, GRB 061007 
decays with a = -1.63 ± 0.11 and GRB 070529 decays with 
a = —1.64 ± 0.14. The other lightcurves in the sample are 
shallower than the reverse shock prediction with > 2a con- 
fidence. The second case arises if i'm,r > Vopt, then imme- 
diately after the peak there is an intermediate stage where 
Q ~ — 0.5, which is followed by q = (3p -I- l)/4. In the 
sample, there are 7 GRBs that before 500s have temporal 
indices consistent with a = —0.5 at 2a confidence. However, 
of these GRBs only GRB 060223a has a decay after 500s 
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Figure 5. The luminosity lightcurves for the 21 lightcurves in the sample with spectroscopic or photometric redshifts. Panel A) shows 
the luminosity lightcurves before correction for the host extinction whereas Panel B) shows the luminosity lightcurves after correction 
for the host extinction. The grey lightcurves are the lightcurves with photometric redshifts. 



(a = — L53 ± 0.60) that is consistent at the 2a level with 
the slowest expected reverse shock decay index a = —1.75. 

The reverse shock in a wind medium is expected to pro- 
duce a muc h steeper decay immediately after the peak with 
a ~ -3.5 (Koba vashi et al.l l2004.'l. Only GRB 050726, has a 
value of Q<5oos, which is consistent to within 2a confidence. 
All other GRBs are inconsistent at > 5(j. 

The inconsistency of the temporal indices of the GRBs 
in this sample to the temporal decay expected during a re- 



verse shock for both types of medium implies that the re- 
verse shock is not the main contributor to the optical emis- 
sion at early times, and therefore is not responsible for the 
rise. Still the reverse shock is expected to occur for all rela- 
tivistic outflows that interact with the external medium. For 
a number of GRBs, the reverse shock may not be observed 
as it can be suppressed by hig h levels of magnetisat ion in 
the outflow (|Zhang et al.ll200d : iGiannios etahl 120081 ') or if 
the forward and reverse shock have comparable energy, the 
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(d) (e) (f) 

Figure 6. Luminosity Distributions. The luminosities in panels (a) to (c) show the luminosities at rest frame 100s, 1000s and 10 ks 
that have not been corrected for host extinction. Panels (d) to (f) show the luminosity distributions at the same epochs, but with the 
lightcurves corrected for host extinction. 



sharp decay in the reverse shock may be masked by the 
flux produced by th e forward shock (|McMahon et al.]|2006l : 
iMundell et alll2007l ). 

5.1.3 Dust Destruction 

If there are h igh levels of exti nction at the beginning of 
the afterglow (|Klotz et al.ll2008l . and references within), the 
hghtcurve produced would be dim and reddened at the be- 
ginning. As the dust is destroyed by the radiation, a chro- 
matic peaJt would be observed as the afterglow brightens and 
becomes less red. The bluer filters would be expected to rise 
more steeply when compared to the red filters as would the 
white filter, because the sensitivity of this filter is skewed to 
the blue. However, the amount of dust destroyed is highly 
dependent on the environment of the GRB, in particular 
to the density and the size of the region surrounding the 
progenitor, and simulations suggest that most of the dust 
destruction is expected t o occur within the first few tens of 
seconds after the trigger (jPerna fc LazzatillioO^ ). Therefore, 
it is unlikely that the UVOT is observing the afterglow while 
dust destruction is occurring. 

However, as the duration of dust destruction and the 
quantity of dust destroyed are only theoretical predictions, 
we must rule out dust destruction using observations. There- 
fore, the 6 UVOT lightcurves with a rise were examined to 
see if the bluer filters, including white, rise more steeply 
when compared to the v filter. The UVOT observed five of 
the lightcurves in white and v (the reddest UVOT filter) dur- 
ing the rise, see Fig.[9l GRB 060607a has the only lightcurve 
where there appears to be a significant excess in v compared 
with the white filter. However, the H band lightcurve given 



in lMolinari et al] |20o3) rises at the same rate as the UVOT 
lightcurve, which was observing in white during the rise. If 
dust destruction was the cause of the rise, the H band would 
be expected to rise less steeply than the UVOT lightcurve. 
The 6th rising GRB, GRB 050730, was observed with the v 
and b filters during the rise and peak. If the peak in this case 
were due to dust destruction, an excess in v compared to b 
would be observed. However, the lightcurves of the v and b 
filters are consistent within la errors. Therefore, there is no 
evidence to suggest that dust destruction is the cause of the 
rise for this GRB or for any of the GRBs in this sample. 

5.1.4 Start of the forward shock 

At the start of the forward shock, a rise will be observed in 
the lightcurves as t he jet ploughs into the external medium 
(|Sari fc Piran|[l999h . The lightcurves for an observer viewing 
within a uniform jet, or within a cone of uniform energy den- 
sity will be the same as thos e observed within an isotropic 
outfiow iGranot et al.l |2002| . and references within). The 
temporal index of the rise will vary according to the thick- 
ness of the shell and the density of the external medium. As- 
suming the synchrotron self absorption frequency Va < Vopt, 
then the thickness of the shell and the density profile of the 
external medium affect the rate at which the lightcurve rises. 
For the thick shell case, the temporal index is q = 1 for a 
constant density medium, or a = 1/9 for a wind medium. 
For the thin shell case, in a constant density medium the 
temporal index is either a — 2 for < i^opt, or a = 3 for 
i^c > Vopt . Lastly, the tempor al index is a = 0.5 if the shell is 
thin and in a wind medium (jPanaitescu fc Vestrand| [2008'l. 
Given the peak time, the Lorentz factor F of the shell 
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Figure 7. Distribution of the rest frame peak times. The red 

filled area represents the GRBs with known rest frame peak times, 
whereas the unfilled area contains the GRBs with known rest 
frame peak times and those GRBs with only upper limits to their 
peak time. Only 21 GRBs, for which luminosity lightcurves were 
produced, are included in this figure. 



tation that the initial Lorentz factor of the jet Fp, wher e 
To ~ 2V{tpeak) (iPanaitescu fc Kuma,ijl2000l : lMeszaroll2006l ). 
of GRBs is > 100 (|Fenimore et al.lll99a ). For the GRBs 
where only an upper limit to their peak time is known, the 
mean value for V{tpcak) is a lower limit, {T{tpeak)) > 230. 
This suggests that the GRBs with observed peaks typically 
have lower Lorentz factors than the GRBs with upper limits 
to their peak times. 

Using the derived values of r(tpeafc), it is possible to 
deduce two more quantities: the isotropic-equivalent mass 
of the baryonic load Mfb = i?/(FoC'^), an d the decelera- 



tion radius Rdec ^ 2ctpeakT{tpeak)^ / + z) l|Molinari et all 
l2007t ). These quantities were determined for each GRB and 
the results are given in Table |6l The mean mass of the 
baryonic load and the mean deceleration radius for the 
GRBs with an observed rise are (M/i,) = 3.8 X 1O-^M0 
and {Rdec) = 2.1 x lO^'^cm and for the GRBs without an 
observed rise the quantities are {Mft) < 1.1 x 1O~^M0 
and (Rdec) < 1-2 x lO^'^cm. The deceleration radii are in 
agreement with Rdec ~ lO^'^cm as predicted by theory 
(iRe es fc Meszaroslll992l ). Therefore the forward shock is con- 
sistent with our observations. 
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Figure 8. Temporal index determined from the lightcurves in 
the rest frame from 150s onwards, shown against luminosity of 
the lightcurves at 100s. 



at the moment of the peak for a constant d ensity medium, 
can be derived using the following equation |Molinari et al.l 
120071 : ISari|[l997l ): 



3E{l + z) 



32nnmpC^ri{tpeak/100s)^ 



1/8 



(3) 



where tpeak is the peak time, 77 is the radiative efficiency 
and n is the density of the medium. Here we will assume 
■q = 0.2 and n = 1 cm~^. However, changing the values of 
rj and n has a minor effect on the final values of r{tpeak) as 
the dependence of T{tpeak) on these parameters is small: 
r(tpeafc) oc {rin)~^^^. For each GRB, the k-corrected en- 
ergy, given in Table |3] was used in the equation above to 
determine r{tpeak) and the resulting values of T{tpeak) for 
the individual GRBs can be found in Table [3] The mean of 
the T{tpeak) for the GRBs with a peak in their lightcurve 
is {T{tpeak)) ~ 180, which is consistent with the expec- 



5.1.5 Off-axis and structured outflows 

A rise may be produced in the lightcurve if the observer's 
viewing angle is Oobs > 0, where Oobs is the observers view- 
ing angle and 9 is the half opening angle of the jet. In the 
case of a uniform jet, the ejecta are released into a cone 
of angle 9 and due to relativistic effects, the emission in 
the jet is beamed as F^^. If 9obs > 9, then the emission 
is strongly beamed away from the observer, but as F de- 
creases, the emission entering the line of sight increases and 
the observed lightcurve w ill rise. The lightcu rve will peak 
when F ~ {9obs -9)'^ (jGranot et al.l I2OO2I '). Larger ob- 
serving angles will view a later peak and the peak magni- 
tude will be lower (|Granot et al.|[2005l ). A structured jet, in 
which the energy per solid angle decreases around a uni- 
form core of angular size 9c, viewed off-axis {9obs > 9c) 
can produce a rise in the optical afterglow, where the be- 
haviour of the rise may vary depending on the distribution 
of energy around the core. In this case, the more diffuse 
the energy per solid angle the slower the rise and the later 
the p eak time (|Panaitescu fc Vestrandl [ioOSi : iGranot et al.l 
|2002| ). The peak time of the lightcurve is also dependent on 
the viewing angle of the jet and the peak will occur when 



5.1.6 Two Component Outflows 

A two component outflow consists of a narrow jet sur- 
rounded by a wide jet. The narrow component will be de- 
noted by a subscript n and the wide component will be 
denoted by a subscript w. Both components move at rel- 
ativistic speeds, but the narrow jet will have a larger F (i.e 
Fu, < Fn) and the wide component will have a larger half 
opening angle than the inner narrow jet (i.e 9^ < 9„). The 
optical emission is ex pected to be produced e ither within 
the w ide component (|de Pasguale et al] l2008l : lOates et al.l 
l2007f ) or mor e traditionally from both the narrow and wide 
components (|Peng et al.ll2005l : iHuang et al.ll2004 ). 
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Figure 9. The early normalized lightcurves for the 6 GRBs with observed rises. GRB 050730 is the only GRB observed to rise that was 
not viewed with the white filter. Key: black star = white, red filled triangle = v, green circle = b, blue diamond = u, light blue cross = 
uvwl, pink hexagon = uvm2, purple square =uvw2 



In the case where the optical emission is only produced 
by the wide component, the rise observed will occur as in 
Section 5.1.4, provided that t he jet is viewe d with in 1.^9^, 
as has been demonstrated by iGranot et al.| (|2002h . At an- 
gles larger than l.bO^, the peak of the optical lightcurve 
will occur when F ~ {9obs — Ow)~ ■ This model of the two 
component outflow could produce the rises observed in this 
sample. 

In the case where the optical emission is produced in 
both the narrow and the wide components, the dominance 
of emission from one component over the other, depends 



mostly on the energy within each component and on the 
vie wing angle of the o bserver (for a more detailed description 
see iPeng et al.ll|2005M . However, in this case as the emission 
is produced in both components, it is likely that on-axis and 
off-ajcis viewers will observe two peaks. This effect is not seen 
in the UVOT lightcurves in this paper and therefore, the 
jet is unlikely to have two components where both produce 
optical/UV emission. 
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5.2 The Late UVOT afterglow 

A correlation has been observed between the observed mag- 
nitude at 400s and tlie decay after 500s. However, there is no 
significant evidence from a Spearman Rank test performed 
between a>i50s, rest and the luminosity in the rest frame at 
100s, for a similar correlation in the rest frame. For an off- 
axis jet, a correlation is expected between the luminosity 
and the decay of the lightcurve when the viewing angle is 
changed, with fainter, shallower afterglows having a larger 
viewing angle jPanaitescu fc Vest rand .2008il . As this data 
do not show a strong correlation of this type, this suggests 
that the outflows are viewed within the half-opening angle 9 
or within a core of uniform energy density 9^. This supports 
the idea that the start of the forward shock produces the 
rises observed in the UVOT lightcurves. However, this does 
not give an indication of the geometry of the jet which may 
be uniform, comprise of two components, or be structured. 



5.3 UVOT afterglow luminosity 

The luminosity distribution found in this work shows no 
evidence for bimodality , which is in contradiction w i th the 
resu lts oflNardini et all l|2006l.l2008l ). lLiang fc Zhangj (|2006l ) 
and iKann et al. ( 2007 ) , who all claim a bimodal distribu- 
tion within their samples. However, the lack of evidence for 
bimodality within t he distributions in t his s ample is consis- 
tent w i th the work of Cenko et al.l l|2008h and lMelandri et al] 
l|2008l ) . ICenko et al.l ( 20081 ) present a luminosity distribution 
at 1000s in the rest frame from a sample of 17 GRBs, with 
known redshift, observed with the Palomar 60 inch tele- 
sco pe. Their distribution shows no evidence for bimodal- 
ity. iMelandri et al.l l|2008l ) produce 3 luminosity distribu- 
tions for 16 optical afterglows observed with the Liverpool 
and Faulkes telescopes. They find a single peaked luminos- 
ity distribution for three rest frame epochs: lOmins, 0.5 
days and 1 day, wh i ch ar e well fit by a log-normal func- 
tion. IMelandri et al.l l|2008l) do not correct their lightcurves 
for host extinction, but as discussed in this paper the cor- 
rection for host extinction appears to have minimal effect 
on the luminosity distributi on. Therefore, i t is possible to 
compare the distributions of lMelandri et aL I diooa^ with the 
distributions produced with this sample. 



5.4 Comparison of the XRT and UVOT canonical 
lightcurves 

There are three segments which are u sually found in X- 
ray lightcurves wit hin the first ~ lO^s (|Zhang et al.l [200^ : 
iNousek et al.|[2006l ). The first segment is a fast, early de- 
cay with —5 < axi < —3, typically ending within 100s- 
1000s after the trigger. The fast decay is thought to be 
caused by the tail of the prompt emission ( Zhang ct al. 2006; 
iNousek et al.l |2006| ) . The second segment is shallow with 
— 1.0 < ax2 < —0.5 (althoug h this range should now be con- 
sidered as —1.0 < ax2 < O.O. lLiang et al.ll2007l ). ceasing be- 
tween 1000s a nd 10000 s and is attribute d to energy injection 
I Zhang et al l 2006; N ousek et al.ll2006l) . The third segment 
decays as —1.5 < axs < — 1 and is expected to occur at 
the e nd of energy injection (jZhang et al ] |2006l : lNousek et al.l 
l2006l ). 

The range in temporal index of the optical lightcurve 



before 500s, taken as the mean plus or minus the dispersion, 
is —1.17 < a < 0.21. This range is clearly inconsistent with 
the first segment of the XRT canonical lightcurve. The range 
of temporal index before 500s is most similar to the range of 
the second segment of the XRT canonical model. However, 
none of the XRT canonical lightcurve segments indicate ris- 
ing behaviour for the X-ray lightcurves. Applying the the- 
oretical interpretations for the individ ual segments of the 
XR T canonical lightcurv e, provided bv lZhang et al.l (|2006l ) 
and iNousek et aTTl|2006h , to the UVOT canonical lightcurve 
suggests that before 500s, the emission producing the opti- 
cal lightcurves is from the forward shock and that a number 
of the UVOT afterglows during this period are energy in- 
jected. The lack of corresponding rising behaviour in the X- 
ray lightcurves, presuming an achromatic rise {um < J^opt), 
suggests that the rise of the forward shock is masked in the 
X-rays, possibly by the contribution of the prompt emission . 
This is consistent with the model of IWillingale et al] (|2007t ) 
who suggest that the steep and shallow decays of the X-ray 
lightcurves are dominated by the prompt emission and af- 
terglow emission respectfully, and that they do not observe 
the rise of the afterglow emission as this is masked by the 
prompt emission. 

In Section 4, the lightcurves after 500s were fitted with 
power-laws and broken power-laws. The GRBs which were 
best fitted by broken power-laws are discussed separately. 
For the optical lightcurves that decay as power-laws, the 
range in decay is —1.22 < a < —0.52. Like the range in de- 
cay before 500s, the range in decay after 500s is most similar 
to the decay range of the second segment of the XRT canoni- 
cal model. Assuming the same reasoning as for the XRT seg- 
ments, this suggests that after 500s the optical lightcurves 
are consistent with emission from the forward shock and 
most of these are energy injected. 

For the four GRBs that are best fit with a broken power- 
law after 500s, the range in temporal decay before the break 
is —0.74 < a < —0.46, which is consistent with the range 
given for the second segment of the canonical XRT model. 
This suggests that before the break the optical lightcurves 
are energy injected. The range in the temporal decay after 
the break is —1.72 < a < —1.34, which is consistent with 
the third decay of the XRT canonical lightcurve, which has 
been suggested to be the decay following the end of energy 
ejection. 



6 CONCLUSIONS 

In this paper we systematically reduced and analyzed a sam- 
ple of 27 GRBs, which met a strict set of selection criteria. 
We note that the temporal behaviour of the optical after- 
glows in the sample is varied, with the greatest variation 
occurring in the early phase of the lightcurves: before 500s 
the lightcurves may rise or decay. The mean for Qf<5oos is 
-OASto ll with a dispersion of 0.69t^;J^. However, after 
500s, all the lightcurves decay. The lightcurves were fitted 
with power-laws and broken power-laws. A broken power- 
law was deemed to be an improvement, if the x^/D.O.F 
decreased and the probability of chance improvement was 
small (< 1%) and in 5 cases a broken power-law was con- 
sidered a better fit. The mean decay index after 500s, when 
including only those that decay with a single power-law, is 
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—0.87^0^(59 with a dispersion of O.SSto'olj- There is a corre- 
lation at 99.8% probabihty, between the magnitude at 400s 
and the temporal decay after 500s, with the brightest optical 
afterglows decaying the fastest. 

We investigated the cause of the rising behaviour in the 
early afterglow and discussed the following physical mech- 
anisms and geometric scenarios: the passage of the syn- 
chrotron frequency Vm, reverse shock, dust destruction, the 
start of the forward shock, the viewing angle of a (possibly 
structured) jet and a two component outflow. The rise in 
the optical lightcurves may be attributed to either the start 
of the forward shock, or to an off-axis viewing angle where 
the observer sees an increasing amount of emission as the 
Lorentz factor of the jet decreases. We also investigated the 
correlation between magnitude and decay after 500s. We de- 
termined that a correlation observed between the magnitude 
at 400s and the decay after 500s is only weakly dependent on 
redshift. A Spearman rank test performed between the lumi- 
nosities at 100s and the decay after 150s in the rest frame did 
not reveal a significant correlation. However, a luminosity- 
decay correlation would be expected for jets viewed off-axis, 
where the more off-axis the jet the fainter and shallower the 
lightcurve. We do not observe a strong correlation of this 
type, suggesting that the optical lightcurves are produced 
by jets viewed on-axis and that the rise observed in the op- 
tical lightcurves is caused by the start of the forward shock. 

We produced luminosity lightcurves for the 21 GRBs in 
the sample with known redshift. The luminosity lightcurves 
were produced at a common wavelength of I6OOA in the rest 
frame. We find that the logarithmic distribution of the lumi- 
nosities at three rest frame epochs: 100s, lOOOs and 10 ks do 
not show evidence for bimodality. Correcting the lightcurves 
for the host extinction increases the mean luminosities of the 
distributions, but does not considerably alter their appear- 
ance and the change in standard deviations of the logarith- 
mic luminosity distributions is no greater than 0.08 for any 
of the three epochs. The lack o f evidence for b i modal ity is 
consistent w i th th e findings of iMelandri et al.l l|2008l ) and 
ICenko et al.1 (I2008D . 

Finally, we compared the temporal behaviour of the 
optical afterglows in this sample with the XRT canonical 
model. We have found that the temporal indices before 500s 
and the temporal indices of the lightcurves after 500s are 
most consistent with the the shallow decaying segme nt of the 
XRT c anonical model. jNousek et al. (2006) and Zh ang et al.l 
l|2006l ) suggest that the shallow decay segment of the XRT 
canonical model is energy injected. This would suggest that 
the optical lightcurves are energy injected as well. The lack 
of rises observed in X-ray afterglows could be due to the 
prompt emission masking them. 
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